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ABSTRACT 


- The  results  of  a  feasibility  study  for  determining  aerosol  size 

distribution  using  the  extinction-sedimentation  Inversion  technique  are 
presented.  It  was  found  that  the  height  of  the  dissemination  chamber, 
may  affect  the  value  of  the  results.  < 


RESUME 


Ce  rapport  presente  les  r£sultats  d'une  etude  de  faisabilitd  de 
1' implantation  d'une  m£thode  optique,  bas£e  sur  la  mesure  du  coeffi¬ 
cient  d'extinctlon,  pour  determiner  la  distribution  granulomitrique 
d'a£rosols  en  sedimentation.  L'etude  a  reveie  que  la  hauteur  de  la 
chambre  de  dissemination  peut  avoir  une  influence  sur  la  valeur  des 
resultats  obtenus. 
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1 .0  INTRODUCTION 


La  dimension  des  particules  est  un  facteur  determinant  de 
l'efficacit6  d'un  aerosol  A  diffuser  la  lumiAre.  Une  ro@me  quantity 
d 'aerosol  diffusera  plus  ou  moins  efficacement  la  lumidre  pour  des 
distributions  granulomStriques  diff6rentes. 

II  existe  plusieurs  m^thodes  pour  mesurer  la  distribution  granu¬ 
lometrique  d'un  aerosol.  II  y  a  d'abord  les  m£thodes  optiques  qui 
regroupent  une  grande  vari£t£  d'instruments  allant  du  microscope  au 
vSlocimStre  laser  A  effet  Doppler.  II  y  a  aussl  les  methodes  non- 
optiques  comme  1 'utilisation  de  grillages  fins,  les  impacteurs  en 
cascade,  1* etude  de  la  sedimentation  de  1' aerosol  dans  le  temps,  et  le 
compteur  Coulter  module  TA  II,  de  Coulter  Electronics  Inc. 

Ce  rapport  presente  les  resultats  d'une  etude  de  faisabilite  de 
1' implantation  d'une  methode  optique,  basee  sur  la  mesure  du  coeffi¬ 
cient  d'extinction,  pour  determiner  la  distribution  granulometrique 
d'aerosols  en  sedimentation  dans  la  grande  chambre  en  forme  de  silo  du 
CRDV. 


La  technique  de  mesure,  appelde  TISE  (Technique  d'Inversion  de 
Sedimentation  et  Extinction),  est  appliquee  A  un  aerosol,  la  poudre  de 
verre  #3419-2  de  Ferro  Industrial  Products  Ltd,  pour  en  determiner  la 
distribution  granulometrique.  Dans  le  but  de  verifier  si  les  condi¬ 
tions  de  validite  de  cette  technique  etaient  respectees,  une  etude  de 
la  masse  d' aerosol  deposee  par  sedimentation  sur  des  papiers-f litres  en 
fonctlon  du  temps  a  ete  faite. 

Les  aspects  theoriques  de  la  technique  TISE  et  de  1' evaluation 
des  sediments  sur  les  papiers-f litres  sont  decrits  aux  chapitres  2  et 
3.  Enfin  les  resultats  obtenus  lors  des  essais  sont  presentes  au 
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chapitre  4.  L'appendice  A  donne  1 'information  relative  a  l'ajustement 
d'une  fonction  math6natique  aux  valeurs  mesur&es  de  la  density  optique. 
L'appendice  B  traite  de  1 'aspect  informatique  de  l'etude. 

Ce  travail  a  effectu€  au  CRDV  entre  mars  et  aout  81  dans  le 
cadre  du  NCP  21B05  "Aerosols". 

2.0  TRANSMISSION  OPTIQUE  AU  TRAVERS  D'UN  AEROSOL  EN  SEDIMENTATION 


L'aspect  th^orique  de  la  technique  de  determination  de  la  dis¬ 
tribution  de  la  grosseur  des  particules  est  dicrit  dans  ce  chapitre. 

On  etudie  en  premier  lieu  la  vitesse  de  sedimentation  des  particules 
sous  1' influence  du  champ  gravitationnel,  puis  on  s’intSresse  A  la 
transmission  d'un  aerosol  monodisperse  en  sedimentation,  pour  en  arri- 
ver  finalement  S  developper  une  expression  analytique  de  la  distribu¬ 
tion  granulometrique  pour  les  aerosols  polydisperses. 

2.1  sedimentation 


Une  particule  spherique  dans  un  milieu  visqueux  non  perturbe 
(absence  de  gradient  de  temperature,  absence  de  vent,  etc.)  tombe  sous 
l'influence  du  champ  gravitationnel  a  une  vitesse  constante  (ref.  1) 
egale  A: 


v 

s 


2ppgr2 

9u 


[1] 


ov)  r  est  le  rayon  de  la  particule  en  cm. 
p  est  la  densite  du  fluide  g/cm3. 

Pp  est  la  densite  de  la  particule  g/cm3. 
g  egale  981  cm/s2. 

U  est  la  viscosite  du  fluide  g/cm.s. 
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Cg  -  1  +  |  (A!  +  A2  e_A32r/A) 


od  l  est  le  libre  parcours  moyen  des  molecules  du  fluide, 

Ai  -  1.257, 

A2  **  0.400, 

A3  *  0.55. 

Le  tableau  I  donne  la  valeur  du  coefficient  C  et  la  vltesse  d 
s§dlmentatlon  de  partlcules  de  diamStre  different. 

TABLEAU  I 

Propri€t£  de  transport  des  aerosols  (r6f.  1) 


Partlcules  sphgriques  dans  1 'air  A  20 °C,  1  atm 


d  (pm) 

P 

C 

v  (cm/s) 
s 

(p  «  1  g/cm3) 

P 

0.1 

2.85 

8.26  x  10~5 

0.2 

1.865 

2  62  x  10”4 

0.5 

1.326 

1.00  x  IO-3 

1.0 

1 .164 

3.52  x  10” 3 

2.0 

1.082 

1.31  x  lO-2 

5.0 

1.032 

7.80  x  10“ 2 

10.0 

1.016 

3.07  x  10"1 

20.0 

1.008 

1.22 

50.0 

1.003 

7.58 

100.0 

1 .00016 

30.3 
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II  apparaft  ainsi  que  pour  les  particules  plus  grosses  que  2  pm 
le  coefficient  Cs  est  A  peu  pris  independant  du  diamitre  de  celles-ci 
et  tend  vers  1.  Par  consequent,  pour  les  particules  plus  grosses  que 
2  ym,  il  est  possible  d'obtenir  une  expression  analytique  simple  pour 
le  diamitre  de  ces  particules  en  fonction  de  leur  vitesse  de  sedimen¬ 
tation,  soit: 


2 .2  Transmission  optique  au  travers  d  ’un  aerosol  monodisperse  en 
sedimentation 


Un  aerosol  monodisperse  est  un  aerosol  dont  les  particules  ont 
toutes  le  mime  diamitre  et  par  consequent,  la  mime  vitesse  de  chute. 

La  figure  1  represente  la  transmission  optique  au  travers  d'un  aerosol 
monodispersi  de  concentration  uniforme.  A  un  temps  tQ,  1'airosol 
est  dissemine  instantanement.  Apris  un  temps  tj  -  t0 ,  il  n'existe 
plus  d'airosol  au-dessus  d'une  hauteur  H-h  (H  et  h  sont  definis  a  la 
fig.  3).  La  vitesse  de  sedimentation  est  donnie  par: 

h 

v  =  - 

S  ti  -  to 

En  substituant  la  valeur  de  vs  dans  l'iq.  2,  on  obtient  le 
rayon  des  particules  constituant  l'aerosol,  soit: 
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Connaissant  le  rapport  Imin/Io»  Imin  et  *o  sont  les 
intensity  lumlneuses  transmises  en  presence  et  en  1 ’absence  de  l'afiro- 
sol,  il  eat  possible  en  utilisant  la  loi  de  Beer  de  connaltre  le  nombre 
de  particules  (par  unitg  de  volume)  qui  ont  £t£  diss£min£es. 

I  ■  I  e"7  ol  T  -  In  I/I  r  6] 

o  o  1  J 


oO  (rgf-  1) 


loo 

T  -  TT  //  K(x  ,  n)  r2n(r)dr  dl  [7] 

00  n 

oO  K^xn»  n)  e8t  l'efficacit£  de  diffusion, 
n  est  1' indice  de  refraction, 
xn  est  £gal  A  2nr/A, 

A  est  la  longueur  d’onde  de  la  lumiere  incidente, 
n(r)  est  la  fonction  de  distribution, 
r  est  le  rayon  d'une  particule, 

L  est  la  longueur  du  chemin  optique. 

La  fonction  K(xn,  n)  peut  itre  eva lu£e  3  l'aide  de  la  thgorie  de  Mie. 
Elle  est  du  type  oscillant  et  tend  vers  2  pour  les  valeurs  £lev£es  de 


Pour  un  aerosol  raonodisperse,  la  fonction  de  distribution  peut 
s'gcrire  2t  l’aide  d'une  fonction  de  Dirac  centrSe  sur  le  rayon  Tp  des 
particules: 


n(r)  -  N5 (r  -  r) 
P 


[8] 


oil  N  est  le  nombre  de  particules  par  unite  de  volume. 
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Apres  integration  et  en  isolant  N,  il  apparaft  que 


Litr2  K(x  ,  n)  2rr2L 
P  n*  p 

pour  les  valeurs  eievges  de  xn. 

II  est  done  possible  de  determiner  la  grosseur  et  le  nombre  de 
particules  d'un  aerosol  monodisperse  en  etudiant  sa  vitesse  de  sedimen¬ 
tation  et  la  transmission  optique  a  travers  celui-ci. 

2.3  Transmission  optique  au  travers  d'un  aerosol  polydisperse  en 
sedimentation 


Une  analyse  du  type  precedent  pourrait  etre  faite  pour  un  aero¬ 
sol  polydisperse  en  utilisant  des  sommations.  Cependant,  Oeepak  et 
Vangham  (ref.  2)  ont  developpe  une  methode  plus  analytique,  relative- 
ment  simple  I  utiliser,  qui  donne  directement  la  fonction  de  distribu¬ 
tion. 


2.3.1  Forme  generale  de  la  fonction  de  distribution 

Considerons  un  aerosol  instantanement  dissemine  et  distribue 
uniformement  au  temps  t  *■  0.  Apr&s  un  temps  t,  toutes  les  particules 
qui  ont  un  diam£tre  plus  grand  que  r(t)  seront  au-dessous  du  niveau  h. 

r(t)  •(?  (Pp  -"<■)  lij  m  (j)‘  l10l 

La  loi  de  la  transmission  optique  (Beer)  devient  une  fonction  du 


temps: 
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I(t)  ,  -T(t) 

~T~  "  e 

o 


[llj 


oO 


r2(t) 

T(t)  -  itL  /  K(x  ,  n)  r2  n(r)  dr  [l2j 

rl  n 


a  un  temps  t  +  At,  on  a: 

r2 (t  +  At) 

T(t  +  At)  ■  ttL  /  K(x  ,  n)  r2  n(r)  dr  [  1 3 J 

ri  n 

En  utllisant  le  th^orSme  de  la  valeur  moyenne  pour  1' Evaluation 
de  ces  deux  IntSgrales  et  la  definition  de  la  d€riv£e,  on  obtient: 

dT(t)  .  T(t  +  At)  -  T(t)  [ 14 J 

dt  At+0  At 

«  xLK(x(t),  n)  n(r(t) )  r2(t)  [l5j 

De  cette  derniire  Equation,  la  fonction  de  distribution  peut 
£tre  isolfie: 

n(r(t))  -  /  *L  K(x(t),  n)  r2(t)  [  16 J 

2.3.2  Equation  de  la  fonction  de  distribution 

Pour  d€finir  compldteraent  1' equation  de  la  fonction  de  distri¬ 
bution  [ 1 6 J ,  on  doit  connaftre  les  derivees  par  rapport  au  temps  de  la 
densite  optique  et  du  rayon  critique  des  particules. 
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Le  choix  d’une  fonction  pour  reprSsenter  la  density  optique 
permet  de  d£finir  analy tiquement  la  fonction  de  distribution.  Pour 
ajuster  les  points  expfirimentaux,  nous  avons  choisi  la  fonction  sui- 
vante: 


T(t) 


T0 

A(eCt  -  1)  +1. 


[17] 


Tq  est  la  valeur  de  la  density  optique  1  t  ■  0,  les  cons- 
tantes  A  et  C  sont  d£termln£es  a  l'aide  de  la  mfithode  des  moindres 
carrSs  aprSs  tranformation  de  l'dq.  17  (voir  appendice  A). 


Les  expressions  pour  les  d£riv£es  sont  obtenues  a  l'aide  des 
6q.  17  et  10. 


dT(t)  -ACT0e 


Ct 


dt  (A(eCt-l)+l)2 


dr(t)  _  1  k  -3/2 

“dt - 2  “  C 


[18] 


L'expression  pour  la  fonction  de  distribution  devient 


n(r(t)  )* 


2ACT0eCt  t5/2 


(A(eCt-l)+l)2  K(x(t),  n)xLa3/2 


[19] 


II  est  pr£f£rable  d'exprimer  la  variable  t  en  fonction  de  r  de 
maniAre  A  Sviter  toute  confusion 


ACToer2  a 

n(r) - - -  [20] 

C— 2 

(A(e  r  -  1)  +  l)2  irLr5 
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Pour  cette  derni£re  Equation,  l'efficacite  de  diffusion  K  a  ete 
flx^e  3  2.  L'gquatlon  20  est  vallde  pour  les  partlcules  de  rayon  plus 
grand  que  1  urn. 

3.0  MASSE  THEORIQUE  DEPOSEE  EN  F0NCTI0N  DU  TEMPS 

La  thfiorle  dScrlte  au  chapltre  2,  p‘ rmettant  de  calculer  la 
fonctlon  de  distribution  granulom£trique  des  partlcules,  suppose  que  la 
sedimentation  des  partlcules  se  fait  dans  des  conditions  ld£ales, 
i.e. : 


-  pas  d'agglomeratlon  par  diffusion  ou  caus£e  par  la  difference  de 
vitesses  de  chute  des  partlcules  de  grosseurs  varlSes, 

-  air  calme, 

-  partlcules  spherlques, 

-  perte  par  diffusion  sur  les  murs  negllgeable. 

Une  faqon  de  detecter  la  presence  de  processus  affectant  la 
vitesse  de  sedimentation  de  1' aerosol,  et  done  de  delimiter  la  region 
de  validite  de  la  technique  etudiee,  consiste  a  mesurer  la  masse  reelle 
deposee  en  fonctlon  du  temps  et  3  la  comparer  3  la  masse  calcuiee. 

Le  calcul  de  la  masse  deposee  par  sedimentation  est  effectue 
pour  une  distribution  discrete,  puis  pour  une  distribution  continue. 

3.1  Calcul  de  la  masse  deposee  par  sedimentation:  distribution 
discrete 


Soit  n  le  nombre  de  partlcules  de  rayon  r  par  unite  de  volume, 
vs  la  vitesse  de  sedimentation,  t  la  variable  de  temps,  h  la  hauteur 


t t ’•£! 
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de  la  colonne  d 'aerosol  au-dessus  du  f litre  au  temps  t  -  0,  et  S  la 
surface  du  papier-filtre. 

Si  le  papier-filtre  est  dficouvert  au  temps  t  -  0,  le  nombre  de 
particules  qui  se  sont  d£pos6es  sur  le  papier  aprds  un  temps  t  est 
donn£  par: 


N  -  nSv  t  H(t  -  h/v  )  [21 1 

p  s  8  1  J 

La  fonction  H  (t  -  h/vg)  est  une  fonction  escalier  (descen- 
dante);  elle  prend  la  valeur  1  pour  les  temps  plus  petits  que  h/v8, 
et  la  valeur  0  pour  les  temps  plus  grands  que  h/v8.  Cette  fonction 
repr£sente  le  fait  qu'aprfcs  le  temps  h/v8  toutes  les  particules  se 
sont  d€pos£es. 

Si  le  papier-filtre  est  dficouvert  3  un  temps  t1 ,  alors 
1' expression  pour  le  nombre  de  particules  d€pos£es  apres  un  temps  t 
est 


Np  -  nSvg  (t  -  t i )  H(t  -  (h/vg  -  ti))  [ 22 J 

Halntenant  si  1 'on  consid£re  des  particules  de  dlam£tre  diffe¬ 
rent,  l'expression  22  devient 

P  fh 

Np  -  (t  -  1 0  SiIi  nlVl  H(t  -  -  tl))  [ 2 3 J 

oO  n^  est  le  nombre  de  particules  par  unite  de  volume,  et  v^  la 
vitesse  de  sedimentation  des  particules  de  type  !• 
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Ce  que  l'on  mesure  sur  le  papier-f litre  n'est  pas  le  nombre  de 
partlcules,  mals  plutbt  la  masse  d£pos£e.  La  masse  d£pos£e  est  §gale 
au  flux  de  masse  a  travers  la  surface  A  multipli£  par  le  temps  £coul£, 
solt: 


m 


8irp  2  .  P 

^(t-tO^f  1  n  r5  H(t  -  -  t  ))  [  24  J 

M  1”1  i  i  Or.  r.r-4  1 


27 


2pp8ri 


2  p  gr2 


ou  la  vltesse  a  iti  remplac£e  par  l'expression  — 2 —  [lj  et  r  est 

le  rayon  des  partlcules. 


3.2  Calcul  de  la  masse  dgposge  par  sedimentation:  distribution 
continue 


Solt  n(r)  une  distribution  granulom£trique  continue,  en  rayon 
des  partlcules,  tel  que: 


N  -  /  n(r)dr  [ 2 5 J 

0 

oO  N  est  le  nombre  total  de  partlcules.  Le  nombre  de  partlcules  d£po- 
s^es  sur  le  papier-f litre  A  un  temps  t  lorsque  celui-ci  a  §t§  d§couvert 
A  un  temps  tj  aprSs  le  d£but  de  la  sedimentation  est  donn£  par: 

00 

Np  -  (t  -  tx)S  f  n(r)  vs(r)  H  (t  ~  ti))  dr  [26j 

La  masse  d^pos^e  sur  le  papier-f litre  est  done: 

8lfp  2  OB  Q 

”  =  -27s-  7T  /  n(r)r5  H  (t  -  tl))dr  [  27  J 

L'  u  0  2ppgr2 

Le  temps  tj  correspond  A  un  rayon  maximum  r^  £gal  A 
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("2p  t  )  »  ce  qui  lmplique  qu'au  temps  ti  les  partlcules  dont  le  rayon 

est  plus  grand  que  r>f  sont  d£jd  d§pos§es  sur  la  plaque  recouvrant  le 
papler-f litre.  Au  temps  tp  od  le  papler-f litre  est  "cachi"  3  nou- 

veau,  toutes  les  partlcules  de  rayon  sup§rleur  3  r  *  (-r - — )  ,  mals 

F  2ppgtF 

lnf§rieur  3  r^,  se  sont  d€pos€es.  Leur  contribution  3  la  masse  de 
sediments  sur  le  filtre  est: 


m(ru>r>rp ^  (V7FJ-^)  n<r>r5dr  t2®] 


rF  8 


oO  v  (r) 
s 


Les  partlcules  plus  petltes  que  rp  se  d&posent  de  fa?on  con¬ 
tinue  dans  l’lntervalle  [  t1 ,  tp]  .  Leur  contribution  3  la  masse 
totale  d€pos€e  sur  le  papier-f litre  est: 

o  rF 

m,  -  p  2  (t  -  tj)  -2-  /  n(r)r5dr  [  29  J 

(r<rp)  27  p  F  u  V  J0  1  1 

La  masse  totale  d€pos€e  sur  le  filtre  s'6crlt  done: 


m  *  m,  N+  m,  .  v  * 
(r<.rF)  (rM>r>rp) 


3.3  Cholx  de  la  fonctlon  de  distribution 


L' Evaluation  de  la  masse  d€pos€e  sur  le  papier-f litre  est  pos¬ 
sible  si  la  fonctlon  de  distribution  granulom€trlque  est  connue.  11 
est  S  noter  que  l'inverse  est  €galement  possible.  La  mesure  de  la 
masse  qui  se  d€pose  en  fonctlon  du  temps  permet  de  determiner  la  fonc- 
tion  de  distribution  (ref.  3). 
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Dans  cette  section,  deux  fonctions  de  distribution  sont  presen¬ 
tees.  La  premi&re  a  §t£  etablie  A  l’aide  d'un  compteur  Coulter,  module 
TA  II  (de  Coulter  Electronics,  Inc)  avant  dissemination  de  la  poudre. 

La  seconde  a  §t£  mesurSe  avec  la  mSthode  TISE. 


3.3.1  Fonction  de  distribution  mesur€e  par  le  compteur  Coulter 


Le  compteur  Coulter  mesure  le  changeraent  de  conductivity  d'un 
electrolyte  place  entre  deux  electrodes  lorsqu'une  particule  passe 
entre  celles-ci.  Le  signal  enregistre  est  proportionnel  au  volume  de 
la  particule.  Apr&s  discrimination  et  manipulation  des  donnees,  l'ap- 
pareil  fournit  les  pourcentages  de  la  masse  de  l'ensemble  des  parti- 
cules  analysees  qui  ont  un  diamitre  plus  petit  que  d. 


Une  mesure  de  la  distribution  initiale  (avant  dissemination)  des 
particules  en  termes  du  pourcentage  de  la  masse  provenant  des  parti- 
cules  dont  le  diametre  est  inferieur  a  d  a  et€  obtenue  pour  la  poudre 
de  verre  #3419-2.  La  density  mesuree  de  ces  billes  de  verre  est  de 
4.2  g/cm3.  La  distribution  granulometrique  en  termes  de  masse  peut 
?tre  representee  par  une  fonction  de  distribution  de  type  Log-Probabi- 

lite  (fig.  2). 


oil  F(m)  est  la  masse  de  l'ensemble  des  particules  ayant  un  diametre  d, 

M  est  la  masse  totale  des  particules, 

o  est  la  deviation  standard  geometrique  par  masse, 
m 

d  est  le  diametre  moyen  geometrique  par  masse. 


50%  des  particules  ont  une  masse  plus  petite  que  celle  d'une  particule 
de  diametre  d,,,. 


II  LA  .IAS SC  AV-tCSSOUS  II  fl 


FIGURE  2  -  Distribution  des  particules  selon  une  echelle  Log- 
Probabilite 

La  courbe  de  la  fig.  2,  telle  que  d£termin£e  par  la  m£thode 
des  moindres  carr€s,  est  celle  qui  repr§sente  le  mieux  les  valeurs 
mesur£es  d  l'aide  du  compteur  Coulter. 

Alnsl,  en  utilisant  1 'expression  mathimatique  reprfisentant  la 
meilleure  courbe  ou  en  proc£dant  graphiquement,  on  trouve  un  diamdtre 
moyen  par  masse  de  10.5  pm.  La  deviation  standard  g£omdtrique  par 
masse  se  calcule  comme  suit  (rgf .  4)  : 

0  _  d84.13Z  _  d  50% 

”  d  50%  d15.87% 

oti  les  diamdtres  indiqugs  correspondent  a  15.87%,  50%  et  84.13%  de  la 
masse.  Celle-ci  est  de  2.42  pour  la  poudre  de  verre  #3419-2. 
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3. 3. 2.1  Distribution  de  type  Log-Probablllt6  en  termes  du 
nombre  de  partlcules 

A  la  section  pr6c6dente  une  distribution  Log-Probabilit£  en 
termes  de  masse  a  et§  ajust§e  aux  mesures  du  compteur  Coulter.  II  est 
n£cessaire  de  transformer  cette  distribution  en  termes  de  masse,  en 
distribution  en  termes  du  nombre  de  partlcules. 

Le  nombre  de  partlcules  de  dlametre  d,  par  unit£  de  volume,  est 
donn£  par: 


F  (d) 


/2n  Log  a 


(Log  d  -  Log  dg)2 
2  Log2o„ 


[33] 


oil  N  est  le  nombre  total  de  partlcules  par  unlt£  de  volume, 

og  est  la  deviation  standard  g€om£trique  par  nombre  de  partlcules, 
dg  est  le  diam^tre  moyen  gfiom£trique. 


50%  des  partlcules  ont  un  dlam^tre  plus  petit  que  dg. 


La  deviation  standard  g6om£trlque  par  masse  est  la  m@me  que  la 
deviation  standard  par  nombre  (r£f.  A).  Le  dlamStre  moyen  g6om£trique 
par  nombre  est  reli€  au  diamStre  moyen  g€om€trlque  par  masse  par  la 
relation: 

Log  dg  *  Log  d^  -  6.9078  Log  [  34] 

En  utilisant  les  valeurs  pour  d^  et  og  de  la  section  pr6c6- 
dente,  on  obtlent  une  valeur  de  1.82  pm  pour  le  diamdtre  moyen  g£om£- 
trique  par  nombre. 
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3. 3. 1.3  Relation  entre  la  fonctlon  F(d)  et  la  fonctlon  n(r) 


La  fonctlon  F(d)  ne  peut  £tre  utllls^e  dlrectement  pour  ^valuer 
lea  dfipSts  sur  les  paplers-flltres.  En  fait,  en  plus  du  changement  de 
variable  d-*-  r,  rappelons  qu'une  distribution  Log-Probablllte  est  le 
r€sultat  d'un  autre  changement  de  variable  (ref.  5): 


done 


n(r)dr 


N _ 

/  2tt  Log  a 


e 


(Log  r-Log  r 

O 

2  Log  o 


)2 

dLog  r 


[35] 


n(r)  -  .  Log  e  [  36] 

Le  nombre  de  particules  par  unite  de  volume  (N)  etant  Inconnu 
jusqu'ici,  11  fera  done  l'objet  de  la  section  sulvante. 

3. 3. 1.4  Density  de  particules  dans  la  chambre  de  dissemination 

L' Evaluation  de  la  masse  d6pos6e  sur  les  paplers-flltres  exlge 
la  connalssance  de  la  fonctlon  n(r).  A  l'aide  du  compteur  Coulter,  on 
a  Itabll  les  proportions  relatives  du  nombre  de  particules,  mals  leur 
density  n'est  pas  connue. 

Une  masse  de  poudre  Mj  est  diss&nin£e  dans  la  chambre  au  temps 
tj)  pendant  une  pgriode  de  30  s.  Lors  de  la  dissemination  et  aprSs, 
pendant  une  autre  p^rlode  de  30  s,  l’intfirieur  de  la  chambre  est  un 
milieu  turbulent  pulsque  les  ventilateurs  de  melange  fonctlonnent. 

Lors  du  melange,  des  pertes  de  particules  survlennent.  II  s'ensult  que 
la  masse  de  poudre  dlss€mln§e  n'est  plus  mals  plut&t  Mp.  La 
determination  de  la  masse  Mp  se  fait  en  mesurant  la  masse  Mfj 
d 'aerosol  d^posee  par  sedimentation  sur  un  papler-f litre  decouvert  au 
debut  de  la  sedimentation. 
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M 

V 


h 

H-h 


[35] 


oil  Sc  et  S  sont  les  aires  du  plancher  de  la  chambre  et  du  papier 
flltre,  les  hauteurs  h  et  H  sont  dffinies  S  la  fig.  3. 


La  density  des  partlcules  est  donn£e  par 


M  /V 

N  -JL-S 

4TTp 


— I  G(r  )r3dr 


[36] 


oO  Vc  est  le  volume  de  la  chambre  et 


G<r>  “  7S-^Ia  -  7 

g 


(Log  r  -  Log  r  )2 
2  Log2  aE 


•  Log  e 


[37] 


3.3.2  Calcul  de  la  masse  dgposge  en  utlllsant  la  fonctlon  de 
distribution  mesurSe  par  la  mSthode  TISE 


Le  calcul  de  la  masse  d£pos£e  sur  les  paplers-f litres  en  utili 
sant  la  fonctlon  de  distribution  mesurSe  par  la  mfthode  TISE  est  rela 
tlvement  simple.  II  sufflt  de  substltuer  la  fonctlon  de  distribution 
n(r)  Stabile  A  l'®q.  20  aux  £q.  28  et  29. 


II  est  a  noter  qu'll  n'est  pas  n§cessalre,  pour  calculer  la 
masse  d^pos^e,  de  determiner  la  densite  de  partlcules,  celle-ci  est 
Incluse  Impllcltement  dans  l'expression  de  n(r). 


4.0  MESURES  ET  DISCUSSION 


Au  chapltre  2,  une  methode  de  determination  de  la  granulometri 
a  ete  presentee;  alors  qu'au  chapltre  3  revaluation  des  dep&ts  de 
masse  par  sedimentation  a  ete  formuiee  matbematlquement.  Dans  ce  cba 
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pitre,  les  r6sultats  exp6rimentaux  et  th^oriques  sont  pr6sent6s  et 
ensulte  discutes. 

4.1  Montage  et  m6thode  exp6rlmentale 

La  fig.  3  illustre  le  montage  utilise.  En  d6pit  du  fait  que  le 
plafond  du  silo  a  une  forme  semi-h6misph6rique,  il  a  ete  consid6r6  que 
l'a6rosol  etait  distribu^  uniform6ment  dans  un  cylindre  de  11.5  m  de 
hauteur.  Un  laser  He-Ne  (632.8  nm,  15  mW  continu)  de  Spectra  Physics, 
et  un  d€tecteur  pyro61ectrique  de  Laser  Precision,  modele  RP313-1,  dont 
l'ouverture  a  et6  r6duite  3  0.5  cm  de  diametre  sont  situes  3  10  m  du 
plafond  de  la  chambre.  Le  d6tecteur  est  couple  a  un  Power  Ratiometer, 
moddle  RK  3441,  de  Laser  Precision;  cet  instrument  est  raccorde  3  un 
enregistreur  graphique. 

Dix  papiers-f litres  Gelman  ayant  un  diametre  de  4.7  cm  sont 
d€pos€s  dans  des  boltes  cylindriques  de  12  cm  de  haut.  Chacune  de  ces 
dix  boftes  est  couverte  d'une  plaque  d 'aluminium.  Neuf  des  papiers- 
filtres  sont  dgcouverts  successivement  3  diff£rents  temps  tfj  a  l'aide 
de  cordes  relives  aux  plaques  d'aluminiura.  Le  dixidme  papier-f litre 
sert  de  reference  et  est  utilise  pour  determiner  l'erreur  sur  la  mesure 
de  la  masse  d6pos6e. 

La  dissemination  de  l'a6rosol  se  fait  a  l'aide  d'azote  a  une 
pression  de  1.03  MPa.  L'aerosol  est  contenu  dans  une  bouteille  de 
plastique  dont  le  col  est  perform  de  petits  trous.  Celle-ci  est  placee 
3  10  m  du  plancher. 

La  sequence  des  6veneraents  au  cours  d'une  experience  est  la 
suivante: 

-  allumer  les  ventilateurs  de  melange  3  t^, 

-  di8s6miner  l’aerosol  pendant  30  s  3  tp 

-  arrSter  les  ventilateurs  au  temps  ts ,  soit  30  s 
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aprds  la  fin  de  la  dissemination, 

-  decouvrir  les  papiers-f litres  au  temps  tjj, 

-  evacuer  le  silo  au  temps  tp. 

Comme  le  rapport  hauteur  sur  diamdtre  des  boftes  est  €lev€,  les 
dipdts  sur  les  papiers-f litres  ne  sont  pas  perturbSs  lors  de  1* evacua¬ 
tion  du  silo  au  temps  tp. 

Beux  experiences  de  sedimentation  ont  ete  effectuees  avec  la 
poudre  de  verre  3419-2  dans  les  conditions  decrites  3  la  section  1.1. 
Le  tableau  II  presente  les  conditions  d’operation  pour  les  deux  expe¬ 
riences. 


TABLEAU  II 


Masse 

Masse  collectEe 

Dur€e  de  fonctionneaent 

Dur£e  de  1' experience 

diapers  £e 

sur  filtre  1 

des  ventllateurs  de 

aprJs  l’arrSt  des  ventllateurs 

a£lange 

de  melange 

(g) 

(8) 

(*) 

(*) 

El 

1310 

0.0217 

30 

5430 

E2 

995 

0.0112 

90 

5430 

4 . 2  Transmission  optlque  et  fonctlon  de  distribution  etablie  par  la 
methode  TISE 


La  fig.  4  represente  la  transmission  optlque  3  travers  la  poudre 
de  verre  #3419-2  se  sedimentant,  en  fonction  du  temps,  dans  les  quan- 
tites  indiquees  au  tableau  II.  La  fig.  5  represente  la  densite  optlque 
dans  la  base  naturelle,  soit  -In  I/IG.  Les  courbes  en  trait  continu 
sont  les  meilleures  qui  representent  les  ensembles  de  points 


El:  1310  g  dlsatalnta,  ptrlode  de  atlange  aprfca 
dlaadalnatlon:  30  a. 

E2:  993  g  dlaalalnta,  ptrlode  de  atlange  aprta 

dlaadalnatlon:  90  a. 


Fond  Jon*  do  dlatr  ibui  Ion  (cm-4) 
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expgrlmentaux  selon  l'expresslon  mathgmatlque  17.  Les  details  relatifs 
au  calcul  des  coefficients  de  la  meilleure  courbe  sont  donngs  en  appen- 
dlce  A. 


Les  fig*  6  et  7  reprgsentent,  sur  des  gchelles  dlffgrentes,  les 
fonctlons  de  distribution  mesurges  par  la  mgthode  TISE  (gq.  20)  en 
utlllsant  les  expressions  mathgmatiques  des  mellleures  courbes  de  la 
density  optlque.  Les  fonctlons  de  distribution  initlale  des  poudres 
avant  dissemination  dgtermi nges  £  l'alde  d'un  compteur  Coulter  et  cal- 
culSes  A  la  section  3. 3. 1.3  ont  §t§  mlses  en  graphlque  en  trait  con- 
tlnu.  On  remarque,  que  pour  les  valeurs  de  diamdtre  glevges,  les 
courbes  des  distributions  sont  parallSles  et  relativement  prSs  l'une  de 
1 ’autre.  Les  courbes  commencent  S  prendre  une  allure  diffgrente  1 
10  |im.  II  est  A  noter  que  pour  les  diamdtres  plus  petlts  que  4.6  pm, 
les  fonctlons  de  distribution  obtenues  par  la  mgthode  TISE  provlennent 
de  1' extrapolation  des  courbes  de  la  densltg  optlque. 

4.3  Comparalson  des  sediments  expgrlmentaux  et  thgorlques 

Les  fig.  8,  9  et  10  reprgsentent  les  sediments  expgrlmentaux 
mesurSs  et  ceux  gtablls  A  l’alde  des  fonctlons  de  distribution  mesurges 
par  le  compteur  Coulter  et  la  mgthode  TISE.  Les  courbes  des  sediments 
expgrlmentaux  se  crolsent  l'une  et  1' autre;  cependant,  ce  crolsement 
n'est  pas  slgnlflcatlf  pulsqu'il  est  £  l'lntgrleur  de  la  marge  d'er- 
reur.  Les  courbes  des  sediments  calculus  A  l’alde  des  fonctlons  de 
distribution  Stabiles  selon  la  mgthode  TISE  se  crolsent  ggalement;  cecl 
est  dfl  au  fait  que  les  fonctlons  de  distribution  se  crolsent. 

Les  fig.  11  et  12  reprennent  les  rSsultats  prgsentgs  aux  fig.  8, 
9  et  10  en  les  regroupant  par  experience.  Les  sediments  obtenus  lors 
de  1 'experience  El,  et  ceux  gtablls  £  l'alde  de  la  fonction  de 


I 


(o ) 


(O) 
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distribution  mesur6e  avec  le  compteur  Coulter  concordent  relativement 
bien  pour  les  paplers-f litres  d6couverts  au  d6but  de  la  sedimentation. 
L'accord  n'est  pas  aussi  bon  pour  les  r£sultats  de  l'explrlence  E2. 

Les  ventllateurs  de  melange  ont  fonctionne  60  s  de  plus  au  cours  de 
1 'experience  E2,  ce  qui  a  favorls£  les  pertes  sur  les  parols  et  1 'ag¬ 
glomeration  des  particules.  Les  particules  §tant  plus  grosses,  elles 
torabent  plus  rapidement;  par  consequent,  la  quantity  de  sediments  est 
plus  eievee. 

L'ecart  relatif  entre  les  sediments  experimentaux  et  ceux  theo- 
riques  (compteur  Coulter)  augmente  avec  le  temps.  Geci  est  dfl  d 
1 'agglomeration  par  diffusion  et  4  1 'agglomeration  produite  par  la 
difference  de  chute  des  particules  de  grosseurs  vari6es.  Alnsl  le 
rapport  entre  les  sediments  experimentaux  et  ceux  theoriques  est  plus 
grand  que  2  apr£s  1000  s. 

L' agglomeration  par  diffusion  d'un  aerosol  monodlsperse  peut 
etre  representee  math6matlquement  par  1' expression  sulvante: 

-  -ovN2  [38] 

oQ  N  est  la  densite  de  particules,  a  est  la  section  efflcace  de  colli¬ 
sion  et  v,  la  vltesse  brownienne  des  particules.  Le  temps  caract6ris- 
tlque  d 'agglomeration  par  diffusion  est  donne  par: 

T  -  (ovN)-1 


Alnsl  pour  N  *  6.4  x  103  cm”3,  d  *  1.82  ym,  et  5  la  temperature 
de  la  piece,  Tc  est  approximativement  1  x  10^  s,  ce  qui  est  une 
periode  de  temps  tr£s  longue.  Cecl  porte  a  conclure  qu'en  raison  de  la 
distribution  relativement  large  des  particules  et  de  la  hauteur  du 
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silo,  1 ’agglomeration  par  diffusion  est  beaucoup  moins  importante  que 
celle  causae  par  la  difference  de  chute  des  particules  de  grosseur 
differente.  A  cet  effet,  les  calculs  de  Friedlander  (ref.  1)  demon- 
trent  que  le  processus  d' agglomeration  par  diffusion  est  dominant  pour 
les  particules  de  rayon  plus  petit  que  1  pm,  et  que  1' agglomeration 
causee  par  la  difference  de  vitesse  de  chute  domine  pour  les  particules 
de  rayon  plus  gros  que  1  pm.  Notons  que  ces  calculs  sont  bases  sur 
1' interaction  de  particules  de  rayon  de  1  pm  avec  des  particules  dont 
le  rayon  varie  entre  0.1  et  10  pm. 

4.4  Discussion 


A  la  section  precedente,  il  a  ete  montre  que  les  sediments  mesu- 
res  et  ceux  calcuies  (compteur  Coulter)  divergent  A  cause  de  1' agglome¬ 
ration  des  particules. 

L'agglomeration  etant  un  processus  cumulatif,  son  effet  est  done 
moins  important  au  debut  de  1' experience.  C'est  ce  qui  explique  que 
les  resultats  anterieurs  5  1000  s  concordent  relativement  bien  avec  la 
valeur  theorique  des  depbts. 

La  methode  TISE  n'est  pas  valide  s'il  y  a  interaction  entre  les 
particules.  II  s'ensuit  que  les  rSsultats  obtenus  A  l'alde  de  cette 
technique  ne  sont  plus  valides  apres  1000  s.  Apres  ce  temps,  toutes 
les  particules  plus  grosses  que  10  pra  sont  deposees.  par  consequent, 
les  fonctions  de  distribution  obtenues  par  la  methode  TISE  ne  sont  pas 
valides  pour  les  particules  dont  le  diamltre  est  plus  petit  que  10  pm. 

5.0  CONCLUSION 


Les  limites  de  la  technique  TISE  de  determination  de  la  granulo- 
metrie  des  particules  ont  §t&  etablies  en  comparant  les  sediments  theo- 
riques  et  les  sediments  experimentaux  sur  des  papiers-f litres  en 
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fonction  du  temps.  II  en  ressort  que  la  hauteur  considerable  du  silo 
augmente  de  faqon  trSs  appreciable  le  temps  de  sedimentation  des  parti- 
cules  et  que,  par  consequent,  1 'agglomeration  devient  importante.  La 
fonction  de  distribution  obtenue  n'est  pas  valide  pour  les  particules 
plus  petites  que  10  pm. 

II  serait  possible  de  determiner  la  granulometrie  des  particules 
plus  petites  que  10  pm  en  deplaqant  le  laser  et  le  detecteur  I  l'extre- 
mite  superieure  de  la  chambre  de  dissemination.  En  effet,  la  hauteur 
de  la  colonne  d'aerosol  au-dessus  de  la  ligne  de  mesure  serait  moins 
grande  et  1 'agglomeration  causee  par  la  difference  de  vitesse  de  chute 
des  particules  serait  moins  importante.  De  plus,  la  duree  de  l'expe- 
rience  etant  plus  courte,  1 'agglomeration  par  diffusion  serait  moins 
grande. 
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En  modifiant  [ A— 1 ]  on  a: 

A  -  A 

Y+  =  -Y°.-A  -f  -  A(eCt  -  1)  +  1  [A- 2] 

En  dfirlvant  l'expression  de  droite 

Y+  -  ACeCt  [A-3] 

et  en  prenant  le  logarithrae  naturel  de  cette  derniSre  Equation  on 
obtlent: 

Z  H  In  Y+  -  Ct  +  D  oQ  D  -  In  AC  [A-4] 


ce  qui  est  une  Equation  lin£aire 
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AI  .2  Calcul  des  coefficients  pour  un  6cart  quadratlque 
min  lml  nun 

Pour  une  Equation  llnSaire  du  type  2  »  Ct  +  D,  le  calcul  des 
coefficients  pour  un  Scart  quadratlque  minimum  se  fait  comme  suit 
(ref.  6): 

S0Ul  -  SiPo  S2P0  “  8lHl 

C - —  ,  D - r  [A-5] 

8082  “  8I  8082  “  8I 

oQ  sq  est  le  nombre  de  points  expSrlmentaux, 

81  est  £  tj , 
s2  est  £  t^, 

PO  est  £  zlt 
Pi  est  £  t^zj. 

Cependant,  dans  le  cas  qul  nous  Intfiresse  11  est  n^cessalre 
d'apporter  quelques  eclairclssements. 

La  valeur  du  coefficient  Aq  est  donnSe  dlrectement  par  les 
coordonnSes  du  premier  point  experimental  solt  (0,ylg)  l.e.  la 
densit€  optlque  au  temps  t  ■  0.  Lorsque  le  temps  tend  vers  l'infini, 
toutes  les  partlcules  sont  dfiposSes  et,  par  consequent,  le  coefficient 
Af  egale  zero, 

Les  €q.  A-3  et  A-4  suggirent  que  l'on  prenne  la  pente  et  finale- 
ment  le  logarlthme  de  celle-ci  pour  llnearlser  1' equation,  Concr&te- 
ment,  on  a 


z 


1 


In  A 

0 


1/yi+l  '  1/yi 

Ci+1  “  C1 


[A-6] 


'  mh  i  ■  ... 
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oO  y^  et  +  i  80nt  les  points  exp£rlmentaux  pris  &  et 
+  11  est  A  noter  que  puisqu'il  a  §t§  n£cessalre  de 

calculer  la  pente,  le  nombre  de  est  &gal  A  n  -  1,  oil  n  est  le 
nombre  de  points  exp£rimentaux. 
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APPENDICE  B 

Informatisation  de  la  m€thode  TISE 
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The  results  of  s  feasibility  study  for  determining  aerosol  size 
distribution  using  the  extlnctlon-sedlnentstlon  Inversion  technique  an 
presented.  It  was  found  that  the  height  of  the  dlsseolnstlon  ch saber, 
say  affect  the  value  of  the  results. 
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